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Дан анализ распределения по вертикали и горизонтали составляющих давления сыпучего материала. Пока- 
зано наличие области максимального давления в центральной части бункера. На основе анализа условия 
равновесия получено выражение для бокового коэффициента. 
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Ведение. В теории взаимодействия сыпучих материалов со стенками емкости важное место от- 
водится определению величины давления материала на стенки емкости. Известно, сыпучий мате- 
риал обладает и свойствами жидкости, и свойствами твердого тела. Модель идеального сыпучего 
материала предполагает наличие силы сухого трения в местах контакта частиц между собой и 
стенками емкости. От величины давления в горизонтальном и вертикальном направлениях зави- 
сит выбор емкости для хранения. Интересным представляется вопрос определения параметров, 
отвечающих за величину давления. 

Теоретический расчет давления в емкостях. В сыпучем материале вертикальное Р, и гори- 
зонтальное Р, давление связаны соотношением: 

Р,=К-Р, (0 
где К- так называемый боковой коэффициент. Это соотношение используется, в частности, и для 
расчета давления на стены элеваторов. Однако, несмотря на то, что соотношение (1) известно 
давно [1], аналитического выражения для этого коэффициента до сих пор не получено. 

В данной работе делается попытка на ос- 
новании аналитического исследования давления Хх 
сыпучего материала в щелевом бункере получить 
в общем виде выражение для бокового коэффи- 
циента _. 

Рассмотрим бесконечно длинную, вдоль 
оси ОХ, емкость с сужающимися книзу стенами, 
т.е. щелевой бункер, с углом наклона стенки к 
вертикали а (рис. 1). 

Система уравнений для определения дав- 
ления в сыпучем материале имеет вид [2]: 
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ду в. д Рис. 1. Щелевой бункер и оси координат 
Р, =К.Р., 


где п, - коэффициент внутреннего трения; у -— объемная плотность сыпучего материала. 


Проинтегрируем первое уравнение в системе (2) по координате у: от -И» до +И/ 
(см. рис. 1), т.е. перейдем от элементарного объема к элементарному горизонтальному слою толщи- 
ной (Е: 


О (3) 


2 





у- 5-42-21”, =(+ 2). 


163 


Технические науки 








Константа интегрирования С по физическому смыслу будет представлять собой верти- 
кальную компоненту силы, действующей на боковые поверхности элементарного слоя. Рассмот- 
рим контакт элементарного, горизонтального слоя с боковой поверхностью щелевого бункера 
(рис. 2). Как видно из рис. 2, к основному элементарному слою АВСР необходимо еще добавить 
два (по одному на каждую сторону бункера) треугольных 
слоя — ВЕС. Если вместо сыпучего материала в бункере на- 
ходилась бы идеальная жидкость, то вертикальная компо- 
нента силы реакции стенки бункера в точности уравнове- 
сила силу тяжести этих элементарных треугольных слоев. В 
первом приближении распространим этот вывод и на случай 
сыпучего материала. 

Следовательно, из сил, действующих на боковые 
стороны элементарного слоя, осталась лишь вертикальная 
компонента силы трения. Силу трения определим как силу 

Рис. 2. Контакт элементарного слоя сухого трения Е» = Ме М, где ре - коэффициент внешне- 
с боковой поверхностью бункера го трения. Очевидно, что сила реакция стенки № будет оп- 
ределяться давлением, действующим на боковую стенку. В 

этом случае константа интегрирования Симеет следующее выражение: 


С=2.и,.М.соза=2.и,.Р.АЕ-соза =2:и,.Р. 








-с0$@. (4) 
0$“ 


Далее, обобщим соотношение (1): 

Р(В)=Р. -1-(-^)-эт В], (5) 
где Р — давление в сыпучей среде по некоторому произвольному направлению; В - угол меж- 
ду этим направлением и осью О7. 

Множитель перед квадратной скобкой — вертикальное давление. Это естественно, так 
как мы рассматриваем движение сыпучего материала под действием сил гравитации. Следова- 
тельно, вертикальное давление будет основным, ведущим фактором, под действием которого 
формируется все поле давления в сыпучем материале. В этом случае, если направление совпада- 
ет сосью 02, то угол В =0и Р(0) =Р, =Р,. Если направление совпадает с осью ОУ, то угол 
В =п/2 и Р(П/2) =Р, = КР». Очевидно, что этот коэффициент должен зависеть от свойств са- 
мого сыпучего материала. 

Подставив в (4) выражение (5), получим: 


С= 2.и,.Р.&=2.и,.4&-Р..[1-(1-К)-соза]. (6) 
Здесь использовано то обстоятельство, что угол между направлением силы реакции М№ и 
осью О2 равен п/2-а. Сучетом (6) уравнение (3) преобразуется следующим образом: 


ОР. 
7.2 :4:2.17, = (+2). > - 4.2.1, +2. и, .4-Р. -1-(-К)-соза]. 
2 
Приведем это уравнение к стандартному виду: 


ФР. _м.-Р.-1-(-&)-соза] 7-8 








= =0. 7 
02 (+2)-("-=-ва) 1+? и 
Здесь использовано соотношение: 
И, =Й-2.юа. (8) 
Решение уравнения (7) имеет вид: 
А 
В 
В =С-С +5.) -_„-СИ+5.2), (9) 


где для сокращения записи введены следующие обозначения: 
у „_ Ш П-(-Ю)-созо] 
1 + и; 1 5 и; | 





ЬБ=юа В= (10) 
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На поверхности бункера давление должно равняться нулю. Из этого условия получим вы- 
ражение для константы интегрирования: 








В _4 
ОЙ ©. 
А-Ь 
При этом уравнение (9) преобразуем к виду: 
А 
И 1 [= |. (11) 
А-Ь 7 
Из выражения (10) следует, что при и, —>0, А->0. При этом Р. > 7 — ‚ что со- 


1+ д; 


ответствует физическому смыслу исследуемой задачи. 
Рассмотрим теперь предел выражения (11) при стремлении стенки щелевого бункера к 


вертикали (а-› 0, Ь-› 0}: 
1 | ИХ }-[77=) ВИ шо 
А Ь—>0 17/4 


50 6550 А-Ь Т-Ь:.2 й 


А А 
Ь Ь 








Таким образом, мы имеем дело с неопределенностью типа 1”. Для ее раскрытия найдем 
сначала предел логарифма этой неопределенности: 
А 
Ъ 


и о ы и] 
ь>0 т ь>0| Б |774 


Для этого разложим логарифм в ряд по степеням |» . = 3 
2 

то ии (1-5-=] о > (5. еее а: 

50| В Ж)| > Ь тои т 


Отсюда вытекает, что: 
А 
_Р.7\Ь 
т [я : еч;[-4 =) 
Ь>0 Ш Ш 


Следовательно, вертикальное давление (при Ь-› 0’) будет стремиться к пределу: 


Шт Р. ее оч И ]. 
ь—>0 А |774 


Подставляя в это выражение значения В и А, получаем: 


В ехр Е , 
Ь—0 и, К 1+1; Ш 


Данное выражение совпадает с хорошо известным решением для прямоугольного 














: К 
силоса [31. При выводе последнего уравнения учтено, что т А= т 5. Таким образом, дока- 
Ь—>0 + и; 
й 
зана преемственность решений для прямоугольного силоса и щелевого бункера. 


Введем теперь в решение (11) зависимость от координаты у. Аналогично [2] эту зависи- 
мость представим в виде: 
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А 

Р _ ВА -Ь.[=-('-у)- вх 1 Е ки ты , (12) 
А-Ь И 

где у - угол естественного откоса. 


Потребуем теперь, чтобы решение (12) удовлетворяло второму уравнению системы (2). 
Найдя первые производные от (12) по уи 2, и подставив их в (1), получаем: 


р И 


вх 8% 

Таким образом, получено аналитическое выражение бокового коэффициента для сыпуче- 
го материала, находящегося в щелевом бункере. Уравнение (13) совпадает с выражением, полу- 
ченным в [2] для сыпучего материала, находящегося в насыпи. 
Тем самым подтверждается предположение о том, что боковой 
коэффициент зависит от свойств самого сыпучего материала и не 
зависит от вида емкости, в которой он находится. 

Однако решение в форме (12) не удовлетворяет нулево- 
му граничному условию на поверхности бункера. Для выполнения 
этого условия, все пространство щелевого бункера разобьем на 
две области (рис. 3). 


Решение в области Г представим в виде [2]: 
Рис. 3. Разделение у.е 
пространства щелевого р. = .2. (14) 


бункера на области 1+ и? 


(13) 


о! \ 








Решение в области П будем искать в виде: 


А 
й — р = == . —р Е с . —-1 . 
ри _ В -Ь Е - У) 5] | |" ь- [И -) г п 
А-Ь И 1+н7 
На линии 2 =0, у = И/ решение (15) обращается в нуль. На границе областей, опреде- 
ляемой уравнением у=И,; 2=(-И/,)-12х (см. рис. 3), первое слагаемое решения (15) 


обращается в нуль, а второе слагаемое переходит в Р/. Таким образом, выполнено граничное 
условие на поверхности бункера. При м--+0 решение (14) преобразуем к виду: 

















РВ [2 -(Т-у)-в1| +8 (И-у)-5х, 
1+1; 1+1; 
или, преобразовав, получим: Р” = т —- - 2, что соответствует физическому смыслу задачи. 
А; 


Итак, мы видим, что решения (14) и (15) удовлетворяют граничным условиям и предель- 
ному переходу к идеальной жидкости. 

График зависимости давления, рассчитанного по формулам (14), (15), от координаты 2, 
при а = 307, у=диу= И! = 3.5 м, представлен на рис. 4. Как видно из рисунка, давление с рос- 
том 2 растет, достигает экстремума, а затем падает до некоторого минимума. Относительное по- 
ложение максимума давления зависит от параметров бункера. 

Таким образом, для идеального сыпучего материала в щелевом бункере найдено верти- 
кальное давление (15). Горизонтальное давление определяем с использованием третьего уравне- 
ния системы (2), выражений (13) и (15). 

Кроме того, исходя из экспериментально измеряемых параметров сыпучего материала, 
получено выражение (13) для коэффициента (К, связывающего вертикальное и горизонтальное 
давление, так называемого бокового коэффициента [1]. Этот коэффициент определяется только 
параметрами самого сыпучего материала (углом внутреннего трения у и углом естественного от- 
коса Хх). В предельном случае, при стремлении к нулю коэффициента внутреннего трения (иде- 
альная жидкость), этот коэффициент стремится к единице, так как в этом случае к нулю будет 


166 


Вестник ДГТУ. 2011. Т. 11, № 2(53) 








стремиться и угол естественного откоса. При этом идеальный сыпучий материал переходит в иде- 
альную жидкость. 





Р, (Па) 


30000 
по образующей 


—’ 
Е. 
й 


20000 


10000 











Рис. 4. Вертикальное давление в щелевом бункере вдоль осевой 
и образующей боковую поверхность линий при различных углах наклона стенки бункера а: 
у = 800 кг / м3; ф = 167; ф = 207; х = 30°; М/ = 3,5 м 


До настоящего времени для определения численного значения бокового коэффициента 
используют эмпирические или полуэмпирические соотношения. Так, например, при инженерных 
расчетах для вычисления этого коэффициента используют следующую эмпирическую формулу: 


к, (16) 
1-и 

где и- коэффициент поперечной деформации. Однако само измерение этого коэффициента для 
сыпучих материалов вызывает большие трудности. К тому же по этой зависимости трудно устано- 
вить предел, к которому стремится боковой коэффициент при переходе от сыпучего материала к 
идеальной жидкости. 

В работе [4] приводятся опытные значения коэффициента бокового давления, в частно- 
сти, для пшеницы К = 0,3 - 0,6, что соответствует приведенным выше расчетам по формуле (13). 

На практике для расчета давления на стены элеваторов используют следующую эмпири- 
ческую формулу для бокового коэффициента: 


к? 45 -" (17) 


Это более реальное выражение для бокового коэффициента. Из этого выражения, в част- 
ности, следует, что при переходе к идеальной жидкости (при ш-› 0) боковой коэффициент стре- 
мится к единице в соответствие с законом Паскаля. Для зерна рекомендуется принимать коэффи- 
циент К = 0,44. Однако в этом случае из (17) вытекает, что угол внутреннего трения должен быть 
равным ц = 23%, что расходится с экспериментальными данными. Выражение (17) нельзя считать 
универсальным из-за отсутствия в нем важнейшей характеристики сыпучего материала - угла 
естественного откоса. В [3] приводится следующее выражение для бокового коэффициента: 

__ № В 
2+1юф.юВ ' 
где [В - угол укладки частиц сыпучего материала. Эта формула получена на основании пред- 
ставлений о дискретной модели сыпучего материала. При этом сыпучий материал моделируется 


абсолютно упругими шарами одинакового диаметра, имеющими правильную геометрическую упа- 
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ковку. Для пирамидальной укладки зерновой пшеницы угол укладки равен 2 = 43°. В (18) угол 
ф - угол внешнего трения. Таким образом, в выражение для бокового коэффициента не входят 
обычно измеряемые параметры самого сыпучего материала: угол внутреннего трения, угол есте- 
ственного откоса. Свойства самого сыпучего материала здесь представлены только теоретически 
введенным углом Д, который экспериментально не измеряется. Поэтому по (18) невозможно оп- 
ределить предельное значение бокового коэффициента при переходе идеального сыпучего тела в 
идеальную жидкость. При ф = 207 коэффициент К, рассчитанный по (18), получается равным 
К = 0,37, что несколько меньше, используемого на практике значения. 

Выводы. Таким образом, можно считать, что расчеты бокового коэффициента дают хорошее со- 
гласие с экспериментом. Формула (13) отражает физический смысл исследуемого параметра, 
удовлетворяет условиям перехода сыпучего материала к идеальной жидкости. 
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